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Nine consecutive days of free-fall, STD data are ana'^yzed in a 
study of the baroclinic structure of the Florida Current. Mean values 
averaged over the nine days are used to reduce tidal aliasing. 

A southward flow is confirmed during the nine day period. Based 
on this data and previous works, the southward flow appears to be of 
a transient nature. 

Ccm.parison of the directly measured and geos trophica lly computed 
current and transport indicates that the Florida Current is essentially 
in geos trophic balance. 
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I. INTRODUCTION 



There have been several previous attempts to determine the 
validity of the geostrophic approximation in the Florida Current*, 

This thesis is a further examination of the baroclinic structure 
of the Florida Current using free-fall, STD instrument data. For 
a period of nine consecutive days in May-June 1969 direct measurements 
of the transport versus depth of the Florida Current between Miami 
and Bimini were made by Dr<> William S. Richardson (NOVA University) 
using free-fall instruments* At the majority of the stations occupied, 
one probe was equipped with the self-contained STD instrument. 

A comparison is made of the average velocity structure determined 
by differentiating the mean transport versus depth curves and of the 
structure determined by geostrophic calculations based on the mean 
density field. To first order, geostrophic equilibrium, i.e. a balance 
between the Coriolis and pressure gradient forces, holds for most 
large-scale oceanic flows. 

The first comparison of the observed and computed velocity fields 
in the Florida Current was made by Wust (1924) in which he used the 
direct measurements of Pillsbury (1890). Considering that the geostrophic 
calculations were based on a density field determined from three 
independent sources, the agreement between the measured and computed 
velocity fields was surprisingly good. The direct measurements of the 
current field by free-fall instruments (Richardson and Schmitz, 1965) 
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have led to more recent comparisons. Broida (1966) used a quasi- 
synoptic density field determined from hydrographic stations, and his 
computations showed a biaxial structure in the Florida Current while 
the direct m.ea sur ements indicated a single axis. The discrepancy was 
attributed to aliasing of the hydrographic da ta by internal tides. 

A time averaged comparison of data taken during the summer months of 
1965 and 1966 was made by O'Brien (1967). Using a mean T-S correlation, 
the density field was determined from the observed temperature only. 
0*Brien*s analysis confirmed the validity of the geostrophic approx- 
imation in the Florida Current. 

The free-fall, STD instrument provides for the first time the 
simultaneous measurement of the transport of all three of the density 
parameters: salinity, temperature, and pressure. By time averaging 

this synoptic da ta , tidal effects are further reduced. Thus, a rare 
opportunity exists to compare accurately the directly measured and 
the geos trophica lly computed velocity fields. 

With an input of geographical station locations, water depths, 
and the values of observed velocity and density versus depth, the 
output of the CHARSECT computer program (Mooers , 1970) includes the 
following baroclinic parameters : 

(1) the thermal wind ratio, which described the degree to which 
a flow is geostrophic, 

(2) the baroclinic stability parameter, which tests the criti- 
cality of the isopycnal slopes, and 

(3) the Richardson number, which describes the dynamic stability 



of the flow. 
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II o METHODS 
A. DATA ACQUISITION 

The free-fall technique yields volume transport per unit width 
versus depth and surface velocity data. When the free-fall instrument 
is equipped with the STD package, a continuous trace of salinity and 
temperature versus depth is also available. The free-fall technique 
employs weighted instruments that fall (attaining their terminal velocity 
of 2 m/sec within a few meters after release) to a pre-selected depth 
where the ballast weights are released, and the instrument returns to 
the surface under its ow^n buoyancy (attaining their terminal velocity 
again within a few meters of ballast release). The precise recording of 
time and position of release and recovery provide the information 
necessary for measuring the depth-dependent transport and the surface 
velocity. At each station one drop is made to the bottom, and one to 
three drops are made to pre-selected depths. Since it is necessary to 
determine the horizontal deflection of the free-fall instrument, the 
navigational system is the controlling factor. The system used is 
Hifix (Decca Navigational System), where the master station is located 
on the vessel and the two slave stations on the western side of the 
Florida Straits. The range of the system is approximately 250 km with 
a precision of +1 meter on the western side of the Straits and +2 meters 
on the eastern side. A sire 11, high speed vessel is used in conducting 
the measurements. 

The data for this thesis were obtained from observations made over 
a nine-day period from 27 May 1969 through 4 June 1969 by Dr, W. S. 






4 



4 



Richardson (NOVA University). A series of thirteen stations vere 
occupied on a section from Miami to Bimini (Figure 1). Despite the 
short crossing time (7 to 8 hours) for a fransect as compared to a 
typical hydrographic transect (20 to 24 hours), the distorting 
influence of tidal motions remains an important factor. In an effort 
to reduce this tidal aliasing, the station times were varied as 
practicable by starting transects at Miami and Bimini on alternate 
days. Tables I and II describe the free-fall instrument data that 
were used. Those observations that included an STD drop are indicated 
in Table II. 

B. DATA ANALYSIS 

The first step in determining the density field was the digitizing 
at ten meter intervals of salinity and temperature values from each 
free-fall STD trace. At each station the mean values of salinity and 
temperature at these intervals were computed for the nine day period. 
These mean values were used as the input to the standard hydrographic 
computer program.^ The output of this program included sigma-t, 
specific volume anomaly, and dynamic depth. The cross stream sections 
of mean temperature, salinity, and (Figures 2, 3 and 4) are plotted 
with the observed current axis superimposed. 

Mean values of observed transport versus depth, rather than 
daily values, were used to compute observed velocity.^ With the volume 

^Mean value is defined as the nine day average. 

^An observed quantity is either a directly measured quantity (such 
as surface velocity) or a quantity (such as sub-surface velocity) 
directly computed from the free-fall measurement of volume transport. 
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Figure 1. Station locations. 
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TABLE I 



Station Number 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


12 


13 


Distance from 
reference point 
(km) 


10 


15 


20 


25 


30 


35 


45 


55 


65 


70 


75 


80 


83 


Maximum depth 
(m) 


135 


285 


345 


350 


340 


665 


765 


795 


830 


790 


735 


580 


395 


Average number 
of drops per 
station occupation 


1 


3 


4 


4 


4 


4 


4 


4 


4 


4 


4 


3 


3 


Total number of 
station occupations 


8 


9 


9 


8 


9 


9 


9 


9 


8 


8 


9 


6 


5 


Total number of 
observations with 
STD 


5 


6 


6 


6 


7 


7 


7 


6 


5 


5 


4 


3 


4 



Free-fall Instrument Data : Miami- Bimini 

Dates : 27 May 1969 through 4 June 1969 



Reference Point ; Lat. 25°44.5'N 
Long. 80°08.8'W 
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STATION NTIMI.ER 


r, 


i 
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4 


5 


6 


7 


8 


9 


10 


n 

11 1 


1 

1 12 


13 


QATE Of 
TRANS EC 1 






















' 








STP 


STD 


STD 


8td 


STD 


std" ' 


' STO 


STD 




STD 


STD - 






Time oi occupation 


0745 


0812 


0900 


0930 


1020 


1115 


1215 


1255 




1^40 


1530 1 




27 May 


Depth of probes 


122 




103 


104 ’ 


102 


173 


173 


171 




340 


170 1 












176 


163 


175 , 


224 


361 


360 


387 




608 


332 ! 












283 


288 


323 


269 


448 


588 ' 


795 




780 


532 1 














325 


3^5 1 


330 


645 


755 








j 












STD 


STD 


1 


STD 


STD 


STD 


STD 


STD 


STu 




STD 


STD 




Time oi occupation 




1428 


1355 




1320 


1240 


1120 i 


1045 


0940 


0903 


0810 1 


0750 


0655 


.28 May 


Depth of probes 




97 


8b 




135 


239 ; 


! 

180 


177 


169 


1 

110 


! 

176 1 


84 


86 








154 


94 




166 


3^:9 1 


341 


346 


338 


357 


354 


320 


181 








270 


329 




335 


4 37 


423 


563 


4 l4 


644 


512 


578 


395 


i 














650 


758 


785 


818 


782 


721 










STD 


STD 


STD 


STD 


STD 


STD ! 


STD 






STD 


STD 




STD 


1 


Time oi occupation 


0806 


0850 


0900 


0942 


1018 


1048 1 


1150 


1245 


1340 


1442 


1512 


1600 


1630 


[29 May 


Depth oi probes 


120 


88 


82 


88 


91 


173 


176 


173 


172 


165 


191 


99 


94 


1 






183 


167 


173 


164 


393 


364 


373 


329 


385 


380 1 


315 


172 








280 


337 


344 


328 


463 


605 


539 


515 


563 


705 


570 


388 










34 0 




335 


657 


747 


773 




780 








' 




STD 


STD 


STD 


STD 


STD 


STD 


STD 


STD 


STD 


STD 


STD 


STD 


STD 




Time oi occupation 


1348 


1324 


1248 


1218 


1148 


1112 


1018 


0918 


0830 


0748 


0712 


0636 


0612 


|30 May 


Depth oi probes 


118 


88 


92 


95 


95 


181 


177 


231 


233 


177 


170 


96 


98 








183 


171 


225 


175 


372 


351 


357 


344 


354 


385 


327 


177 








270 


265 


330 


326 


494 


598 


520 


525 


523 


695 


570 










1 


320 


350 


340 


636 


748 


780 


806 


776 








































Time oi occupation 


lono" 


0730 


0810 


0930 


1000 


1030 


1115 


1200 


^ 1300 


1330 


1350 


1415 


1430 


51 May 


Depth of probes 


133 


95 


99 


93 


102 


91 


96 


95 


93 


92 


90 


90 


91 








179 


235 


185 


134 


172 


166 


160 


160 


208 


167 


372 


133 








277 


277 


316 


310 


607 


600 


613 


673 


610 


474 


573 


395 


j 








326 


350 


336 


662 


754 


782 


814 


765 


733 
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1150 
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1120 




1040 
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0900 
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0710 
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1 June 


Depth of probes 


94 


278 


97 




90 


94 


98 


95 


103 


95 


86 


93 












160 
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177 
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173 


152 


176 


721 


290 
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763 
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812 
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STD 
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2 June 


Depth of probes 


125 


87 


89 


84 


89 


91 


89 


115 


86 


87 


158 












169 


176 


176 


228 


165 


165 


242 


172 


171 


310 












275 


250 


328 


278 


443 


480 


483 


509 


465 


679 














340 


340 


340 


658 


548 


780 


800 


737 
























750 
























STD 


STD 


STD 


STD 


STD 


STD 


STD 


STD 


. STD 


STD 




Time of occupation 


1712 


1700 


1625 


1555 


1530 


1450 


1400 


1312 


1215 


0840 


0700 


' 0615 


0530 


3 June 


Depth of probes 


128 


92 


89 


99 


84 


93 


145 


96 


93 


96 


81 


240 


34 








181 


182 


171 


178 


243 


167 


161 


187 


144 


169 


288 


: 183 








282 


255 


282 


291 


428 


510 


488 


236 


182 


367 


580 


395 
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330 


645 


750 


788 


466 
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610 


























810 


780 












STD 
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STD 


STD 


STD 
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0842 


0906 
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1012 


1042 


1118 


1236 


1324 


1430 










^ June 


epth of probes 
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96 


92 


85 


154 
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230 


1 228 
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; 363 
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Figure 2. Mean temperature (C^) section 
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Figure 3. Mean salinity (%q) section 
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Figure 4 



Mean sigma-t section 
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transport available at several depths, a transport versus depth curve 
was drawn for each occupation of each station. Values were read from 
each curve at 50 meter intervals; these values were then averaged to 
obtain the mean volume transport curves. By differentiation of these 
curves, the mean velocity versus depth profiles were obtained. 

Using the method of least squares, second, third and fourth order 
polynomials w’ere fitted to these transport curves. Two constraints 
were placed on the polynomials : 

1. The transport was forced to equal zero at zero depth (sea 
surface . 

2. At zero depth, the derivative of the polynomials was forced 
to equal the observed mean surface velocity. 

An examination of the various orders of polynomials showed that the 
third order generally had the best fit and resulted in the most realistic 
profile. The third order polynomial was also used as the best fit in a 
similar treatment of free-fall data (Richardson, Schmitz and Niiler, 
1969). Comparing the velocity versus depth profiles at stations 5 and 
7 for second, third and fourth order polynomials (Figure 5), it appears 
that the third order polynomial yields the most realistic profile. 

The methods described above for determining a mean transport curve, 
and, subsequently, the velocity versus depth profiles, w’ere not the 
only ones attempted. After several tries, it was found that 50 meters 
was the minimum increment that could be used for fitting a polynomial 
to a mean transport curve. 

Another method tried was a least squares polynomial fit (with the 
same tw’o constraints) to the raw transport data, i.e. the actual values 
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Figure 5. Observed velocity profile for stations 5 and 
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of transport available at several depths,, This method resulted in 
less realistic velocity profiles, particularly in stations such as 
11, 12 and 13 where fewer observations w'ere made. 

At a sample station, direct differentiation of the mean transport 
curve resulted in a velocity profile that was comparable to the one 
obtained by differentiation of the polynomial- The method used in 
this thesis cannot be said to be the "best" method, but it does yield 
velocity profiles that are reasonable when compared with previous works. 

The resulting mean axial velocity structure across the Florida 
Straits is shown in Figure 6. The current axis of this nine day mean 
velocity field was superimposed on the temperature, salinity and 
sigma-t sections (Figures 2, 3 and 4). Of particular note is the 
southward flow near the bottom between stations 2 and 8. This southerly 
flow, which is discussed in Chapter IV, is in agreement with Wust (1924), 
bottom photos of ripplemarks by Neumann (1970), and the direct measure- 
ments by Duing and Johnson (1971) and Duing (1971). 

The geos trophic currents w-ere computed according to the classical 
dynamic method. As derived by Sandstrom and He lland-Hansen (1903), 
the relative velocity between two isobaric surfaces, under the 
assumption of geostrophic flovr, can be expressed as 

1 0 P2 P2 

V2-V1 = ^ f / QAdp - / 

Pi Pi 

where 

f is the vertical component of the Coriolis parameter, f = 2o)sin0 

and U) is the angular rotation of the earth and 0 is the latitude; 

Ax is the horizontal distance between stations A and B; 
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Figure 6. Mean observed axial velocity (cm/sec). 
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Oi is the specific volume, = ^35,0,p ^ ^ specific 

volume anomaly; 

and, V\^2 pressure on isobarib surfaces 1 and 2. 

The dynamic depth, D, at isobaric surface p^^, referenced to an isobaric 
surface pg, is 

Pi 

D = / adp (II-2) 

PO 

By substitution of (II-2) into (II«1), the relative velocity between 
isobaric surfaces 1 and 2 becomes: 

V9-Vi (^2 )a’ (^i)a^ “ C(D2)b“(Di)b] (II~3) 

f _ 

With the field of mass known, (II-3) can be used to compute the 
relative current field (Neumann and Pierson, 1966). To determine the 
absolute velocity field, there must be at least one level, the reference 
level, where the absolute velocity is known. The establishment of a 
reference level is usually difficult and speculative. Computation of 
the absolute velocity field is possible if good quality direct measure- 
ments are available. 

The following direct and indirect methods were used to determine 
the reference level across the Florida Straits: 

(1) Determination from the observed data of the depth where the 
observed velocity equals zero. (direct) 

(2) Setting the observed and computed velocities equal at the 
surface. (direct) 

(3) Defant’s AD method. (indirect) 

The depth of no motion is the depth where the mean current component 
normal to the section is zero. The use of direct measurements to 
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determine the depth of no motion is best suited to a region where the 
mean current is strong. The observed velocity versus depth curves 
yields the depth of no motion, if one exists, for each station. Knowing 
the depth of no motion, the absolute velocity field can be computed. 

Alternatively, the directly measured surface current can be used 
to convert the relative velocity field into an absolute one. A ref- 
erence level results by setting the geostrophic surface velocity equal 
to the observed surface velocity. On the anticyclonic side, equating 
the two surface velocities did not indicate a depth of no motion in 
some cases, i.e. northward flow filled the entire water column.^ 

For strongly stratified fluids, one of the most reliable indirect 
methods is that proposed by Defant (1941). The Defant method consists 
of comparing relative dynamic depth differences at given isobaric 
surfaces between two stations and to determine at what level the dynamic 
depth difference is constant. This process defines a level with a 
constant horizontal pressure gradient, which is hypothesized to be the 
depth of no motion. 

The differences in the dynamic depths, computed relative to the 
sea surface, were compared for various combinations of station pairs. 

In the majority of the comparisons there was a definite layer where 
Ad was constant. Near the boundary on the anticyclonic side, several 
of the comparisons did not show a distinct level of constant Ad. 

Figure 7 shows typical Ad profiles for the cyclonic side, the middle, 
and the anticyclonic side of the Florida Current. For these profiles 

^Anticyclonic side in the Florida Current is where the relative 
vorticity is negative, i.e. the eastern side. The cyclonic side is 
where the relative vorticity is positive, i.e. the western side. 
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igure 7. Mean d 3 mamic depth difference (AD). 
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the layers over which Ad is constant are 60, 70 and 5 meters respec- 
tively. 

Using the three methods outlined above, a section showing the 
depth of no motion was computed (Figure 8). The extension of the 
reference level from m.ethod (2), Vgeos trophic = Vobserved at the 
surface, into the eastern boundary occurs where northward flow filled 
the entire water column. The doubtful areas of the reference level 
from method (3), the Defant method ^ are indicated by dots superimposed 
on the solid line. It is clear that the differences between these 
m,ethods are not great. Typically it is a 25 meter difference west of 
station 8 and a 50 meter difference east of station 8. 

Different geostrophic calculations were made using various 
reference depths and various station pairs. In some cases the reference 
depth was below the maximum depth of one of the stations due to sloping 
boundaries. In these cases it was necessary to extend the dynamic 
calculations below the bottom. For these calculations, the dynamic 
depth profile was extrapolated to the necessary depth. 

Downstream velocity sections were computed using a depth of no 
motion as determined by Defant *s method (Figure 9) and by matching of 
observed and computed velocities at the surface (Figure 10). Both 
figures have the observed velocity field superim.posed . Compared to 
the geostrophic isotachs, the observed isotachs^, except near the surface, 
are deeper and skewed to the east. In both figures, the best agreement 
between observed and geostrophic occurs on the cyclonic side. 

Two additional indirect methods are shown in Figures 11 and 12. 

A velocity section was computed using the bottom as the depth of 
no m.otion. and compared with the observed field (Figure 11). The 
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Figure 8. Depth of no motion 
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Figure 9. Mean axial velocity (cm/sec), observed and geostrophic 
(depth of no motion by Defant*s method). 
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Figure 10, Mean axial velocity cm/sec;, observed and geostrophic 
(^geos trophic ~ ^observed surface). 
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Figure 12. Mean axial velocity (cm/sec), observed and geos trophic 
(depth of no motion = 400 meters). 
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biaxial structure is similar to the results of Broida (1966). Broida *s 
biaxial structure was attributed to internal tidal aliasing. In this 
case the structure is due to the use of a non-rea listic depth of no 
motion on the cyclonic side. 

To illustrate the effect of a constant depth of no motion, a 
section was computed using 400 meters as the depth of no motion (Figure 
12). Again the biaxial structure is a result of an unrealistic reference 
depth. 

In each of the geos trophic sections computed, the following com- 
binations of station numbers w^ere used: 1-3, 2-4, 3-5, 4-6, 6-8, 7-9, 

9-11 and 10-12. All isotachs w*ere extrapolated to the boundary at the 
surface or interior. 



C o ERRORS 

The errors associated with the free-fall instrument data have two 
sources: the system errors and the possible interpretative errors 

arising from the assumptions that conditions are stationary over the 
time and space scales of a run. The system error consists of navi- 
gational and timing errors. The Hifix has an error of +1 meter on the 
western side of the Straits and +2 meters on the eastern side. The 
timing error is +1 second. The system error yields an error estimate 
for the surface velocity of +3% and for the transport of +5% (Richardson 
and Schmitz, 1965). For the short run time (1-10 minutes) and the short 
horizontal distance of a run (10-500 meters), experience has demonstrated 
that fluctuations large compared to the mean axial flow do not occur 



over these scales in the Florida Current. 
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A transect is completed in about 8 hours. As stated previously, 
mean values vice daily values were used to reduce the effect of tidal 
aliasing- In addition, the station times were varied as much as 
practicable to reduce tidal aliasing. Examination of the times of 
station occupation over the nine day period (Table II) reveals that 
there is little variation in those times for the stations in the 
middle of the Straits. The diurnal tidal motion may therefore have a 
greater aliasing than the semi-diurnal. Furthermore, the solar 
diurnal may have a greater aliasing effect than the lunar diurnal. 

Using the times of station occupation and values of the various 
tidal harmonic constants from Florida Current data (Smith, Zetler and 
Broida, 1969), the tidal aliasing error in the observed velocities is 
+37o for the axial velocity. The error is approximately 257o for the 
cross-stream velocity. Appendix A gives the procedure used to calculate 
these errors. 

The values of depth, temperature and salinity vrere read directly 
from the STD trace Vvithin the limits : 
depth, +0.5 meters; 
temperature, +0.025^; and 
salinity, +0.0^o • 

The geostrophic velocity error is a sum of the relative errors in 

depth, horizontal separation, and density. The relative error in the 

density of +^0% is more than one order of magnitude greater than the 

relative errors in depth and horizontal separation. Hence, the 

relative geostrophic velocity error is +207o. 

A comparison of O' and p. at the deepest station, number 10, 

^ t '^in-situ ^ 




t 




ST 
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showed that the approximation of the density field by resulted in 
negligible error. 

For an indication of the variability of the data, the standard 
deviations for temperature, salinity and sigma-t are plotted (Figure 13), 
The area of maximum variance in these parameters is in the pycnocline 
on the cyclonic side. 
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III, RESULTS 

A. COMPARISON OF DIRECT MEASUREMENTS 
WITH GEOSTROPHIC MEASUREMENTS 

Other than on the anticyc Ionic side, where the Defanc method is 
in doubt, there is a good agreement between the directly and indirectly 
determined depths on no motion (Figure 6) o From this figure, a ”best‘' 
depth of no motion, called a hybrid depth of no motion, is chosen and 
used for the final comparisons of observed and computed fields. On 
the cyclonic side, the hybrid depth of no motion is the level determined 
by the Defant method. On the anticyclonic side, it is the bottom. 

Since the hybrid depth of no motion nearly coincided with the depth 
where the observed velocity equals zero, the latter could have been 
used as the best depth of no motion. The hybrid depth was chosen for 
the following reasons : 

(1) it does closely approximate the depth of no motion determined 
by the observed velocity field, and 

(2) its use sets a precedence for determining a depth of no 
motion in the Florida Straits when no observed velocity data 
are available. 

The computed surface velocity distribution using both the Defant 
and hybrid reference depths is compared to the observed surface velocity 
distribution in Figure 14. On the anticyclonic side the agreement 
between computed and observed surface velocities is poorer than on the 
cyclonic side. East of station 9, the surface velocities computed with 
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Figure 14. Mean surface velocity (cm/sec), observed and geos trophic (hybrid depth of no motion). 
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the hybrid depth shows less agreement with the observed surface 
velocities than those computed using the seemingly less reliable 
depth found by Defant's method. 

Using the hybrid depth of no motion, a comparison of the computed 
and observed velocity sections is made (Figure 15), As with the 
surface velocities, the best agreement is on the cyclonic side of the 
stream. Compared to the geos trophic velocity, the observed isotachs 
are deeper, except near the surface, and the observed current axis is 
skewed to the east. The agreement between the position of the observed 
current axis and the position of the geostrophic current axis improves 
as depth increases. 

The geostrophic, using the hybrid depth of no motion, and observed 
velocity profiles are compared in Figure 16 for each station where 
both computed and observed values were available. The geostrophic 
subsurface maximum at station 11 is similar to the observed subsurface 
maximum by Duing and Johnson (1971), The closest fit of the absolute 
values of the curves occurs on the cyclonic side. However, comparison 
of the vertical shear shows the best agreement on the anticyclonic side. 

Because of the uncertainties in determining the observed velocity 
fields by differentiation of the mean transport curves, an informative, 
if not more accurate, test of the validity of the geostrophic approx- 
imation is the comparison of transport curves. Figure 17 shows 
transport per unit width curves, observed and computed. A representative 
station from the cyclonic side, the middle and the anticyclonic side is 
used for the comparison. The previously established pattern of close 
agreement of absolute values on the cyclonic side and poorer agreement 
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STATION . NUMBER 12 3 4 5 6 



9 10 II 12 13 




Figure 15. Mean axial velocity (cm/sec), observed and geostrophic 
(hybrid depth of no motion). 
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on the anticyclonic side is clearly illustrated. 

The thermal wind relation derived from the assumption of geo- 
s trophic and hydros ta.tic equilibrium is 
f?2 = 

Po 

where 

V 2 is the derivative with respect to depth of the mean axial 

velocity, 

is the derivative with respect to horizontal distance of the 

mean density, 

pQ is the reference density, 
and, g is the gravitational acceleration. 

Using the observed velocity and density data, the ratio of left 
hand side to the right hand side describes whether or not the thermal 
wind relation is satisfied. If the ratio is unity, the thermal wind 
relation is satisfied. If the ratio is greater (less) than one, the 
observed shear is greater (less) than the geostrophic shear. 

From Figure 18, the areas of observed southward flow and the area 
east of station 9 are where this ratio differs the greatest from unity. 

B. FEATURES OF THE FLORIDA CURRENT 

To extract additional information about the baroclinic structure 
of the Florida Current, the free-fall data were analyzed by studying 
the following: 

(1) Baroclinic stability parameter 

(2) Richardson number 



(3) Mean T-S curves 



37 



STATION NUMBER I 
DISTANCE (km) 10 



2 3 4 

20 



5 6 7 

30 40 



8 

50 



9 

60 



10 II 12 13 

70 80 




Figure 18. 



Thermal wind ratio section. 



(4) Net transport between selected isopycnals 

(5) Total downstream transport 
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(6) Cross “Stream flow 

The results of these studies are discussed below. 

(1) The baroclinic stability parameter is the ratio of the 
slopes of an isopycnal, S, to its critical value, Sc (Mooers , 1971). 

When S > Sc, baroclinic instability, or hydrodynamic disequilibrium, 

can occur. The baroclinic stability is relatively low in the pycnocline 
on the cyclonic side (Figure 19). 

(2) The gradient Richardson number is expressed as 

Ri = — . 

(^z) 



where 

N^ = ^g is the Vaisa la-Brun t frequency. 

Ri<l implies dynamic instability for the flow, and Ri>l implies 
stability. From Figure 20, the area where the dynamic stability is 
the lowest is near the bottom on the cyclonic side. The dynamic 
stability is the greatest near the surface on the anticyclonic side. 

(3) The mean T-S diagrams for stations 2, 6 and 10 (Figure 21) are 
in agreement with the T-S diagrams presented by Wennekens (1959) for 
the Florida Current. The cross channel distribution of water mass 
properties described by Wennekens is illustrated in the T«S curves. 
Station 2 is representative of the Continental Edge Water, station 6 
is representative of the Transition Zone Water, and station 10 is 
representative of the Yucatan Water. The Yucatan, or Caribbean, Water 
is identified by its well defined salinity maximum. The great reduction 
in the intensity of this salinity maximum is the conspicuous feature of 
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Figure 19. 



Baroclinic stability parameter section. 




Figure 20. Richardson number section (isolines have non-uniform 
spacing) . 
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T-S curve, stations 2, 6 and 10. 



Figure 21. 
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the Edge Water. 

The water that flow's through the Florida Straits originally comes, 
in large part, from the southern half of the North Equatorial Current 
and from a branch of the South Equatorial Current. This w’ater flow’s 
through the Caribbean, and then, w’ithout mixing with the w’aters 
endemic to the Gulf of Mexico, passes through the Florida Straits in 
very nearly its original state. Because this water has acquired a 
large admixture of Antarctic Intermediate Water at mid-depths from 
the South Atlantic, there is a salinity minimum (betw’een 600 and 800 m 
depth) in the w’ater exiting out of the Florida Straits (Stommel, 1966). 
This salinity minimum is present in stations 6 and 10. 

(4) The net axial transport across the Straits and between 
isopycnals was determined in the vicinity of the southw’ard flow’, A 
net transport southw’ard would favor the existence of a southward 
undercurrent rather than a large scale eddy. In all combinations of 
isopycnals, the net transport was alw’ays greater to the north. This 
result admits the possibility of a large scale eddy being the cause 
of the southw’ard flow. An example of the results are shown below. 



Isopycnal 


Interval 


Indus ive 


Stations 


Net Transport 


27.2 - 


27.3 


2 - 


9 


7.5 X 10^ m^/sec North 


27.4 - 


27.5 


6 - 


9 


3 X 10^ m^/sec North 



(5) The net downstream transport value at each station was 
integrated across the channel (Figure 22), The total mean volume 
transport of 33.4 x 10^ m /sec compares favorably with the value 
obtained by Richardson and Schmitz (1968) of 32.2 x 10^ m^/sec for data 
averaged over the period of May-June 1965. 



STATION NUMBER I 
DISTANCE (km) 10 
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Figure 22. Mean observed downstream transport (m^/sec) and surface velocity (cm/sec) with standard 
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The volume transport for each of the nine days was computed 

separately and is shown in the table below. The average of these 

b o 

nine values is 33.4 x 10^ m'^/sec. 



Da te 


Volume transport 
(x 10^ m^/sec) 


27 May 


30.3 


28 May 


33.3 


29 May 


34.7 


30 May 


33.2 


31 May 


36.8 


1 June 


31.2 


2 June 


36.5 


3 June 


31.3 


4 June 


32.9 



Continuous electrode potential measurements of transport based on 
two electrodes located at stations 2 and 3 showed good agreement with 
the free-fall measurements during the nine day period. The directly 
measured and electrode measured transport differed by less than 107o, 
but tidal aliasing of the free-fall measurements precludes any firm 
conclusion (DeFerrari, 1970). 

(6) The cross-stream data is not of sufficient quality for a 
detailed analysis of the mean cross-stream velocity structure. A 
general pattern can be ascertained by the differentiation of the mean 
transport curves. On the cyclonic side, the velocity is westward except 
for a mid-depth layer (50-70 meters thick) of eastward flow. On the 
anticyc Ionic side, the velocity is eastward in the upper 400-500 meters 
and westward below this depth. 



The mean cross-stream transport is shown in Figure 23, where 
positive values indicate eastward transport and negative values, 
westward transport. There is westward transport in the westward side 
of the Straits, eastward flow in the eastern side of the Straits, and, 
hence, an area of divergence near the center. The mean patterns are 
statistically significant at the 95% level. 



STATION No. I 2 ^ 4 5 6 7 8 9 10 il i2 13 

DISTA NCE ( k,t.) 1,0 . 2.0 , 3.0 , 4.0 , 5,0 , 6,0 , 7,0 , 9, 0 . 
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Figure 23. Mean observed cross-stream transport (m^/sec) with standard errors 
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IV, SUMMARY AND CONCLUSIONS 
A, SUMMARY 

Nine consecutive days of free-fall, STD data, taken by Richardson 
during late May and early June of 1969, were analyzed in a study of the 
baroclinic structure of the Florida Current, The directly measured 
transport was used to obtain the mean velocity field, and the directly 
measured temperature and salinity profiles were used to obtain the mean 
density field. Mean (ioe. averaged over the nine day period) vice 
daily values were used to reduce the tidal aliasing problem. 

The mean axial velocity vras southward beneath the Florida Current 
during the nine day period. This indication of southward flow was 
consistent w’ith Hurley (1963), Neumann (1970), Duing and Johnson (1971), 
and Duing (1971). Richardson, Schmitz and Niiler (1969) did not find 
this southward flow as permanent feature when more extensive data were 
analyzed. 

For the geos trophic calculations, three methods were used to 
determine the depth of no motion: the depth where the observed 

velocity equals zero, equating the geostrophic and observed velocities 
at the surface, and Defant's method. On the cyclonic side, there was 
good agreement between the depths of no motion determined by the three 
methods, and on the anticyclonic side the agreement was not so good. 
Likevrise, on the cyclonic side Defant's method showed a definite depth 
of no motion, but on the anticyclonic side the results vrere less 
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conclusive. 

Using the hybrid depth of no motion, the geos trophic velocity 
field was computed. The comparison of the observed and computed 
velocity fields showed good agreement on the cyclonic side and much 
less agreement on the anticyclonic side. The same pattern of agreement 
between absolute values existed in the comparison of surface velocities 
and of velocity and transport profiles at several stations. The thermal 
wind equation was shown to be a good approximation, i.e. the flow was 
essentially geostrophic throughout the Straits. The greatest dis- 
crepancies occurred in the area of the southward flow and in the 
easternmost part of the anticyclonic zone of the Straits. The dis- 
crepancy may be related to the limitations of the observational and 
data analysis techniques. 

Also, the study of the standard deviation of temperature, salinity 
and showed a zone of high variability on the cyclonic side. 

The study of stability parameters showed an area of low dynamic 
and baroclinic stability on the cyclonic side. 

B. CONCLUSIONS 

(1) The validity of the free-fall, STD measurements was established 
by comparison with previous total transport values. 

(2) The Defant method was a valid method for determination of a 
depth of no motion on the cyclonic side. The bottom was a reasonable 
depth of no motion on the anticyclonic side. Thus, absolute geostrophic 
velocities can be computed independent of a direct method for depth 

of no motion determination. 

(3) The Florida Current was essentially in geostrophic balance. 
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(4) A southerly flow beneath the Florida Current was confirmed 
during the nine day period. The T-S curve of the mid-channel southward 
flow has shown the salinity minimum that is characteristic of the 
Antarctic Intermediate Water. If a steady countercurrent existed, 
different T-S curves would be likely. In the vicinity of this 
countercurrent, the net transport between isopycnals was substantially 
to the north. The lack of net southward transport, and of a distinction 
in the T-S correlations, precluded resolving whether a large eddy or 

a steady countercurrent exists on the basis of the present data. 
Richardson, Schmitz and Niiler (1968) showed that, over a longer time 
span, the north component of velocity fills the whole channel, which 
implies a transient nature for the southerly flow. The velocity 
profile analysis by Duing and Johnson (1971) and Duing (1971) gave 
definite indications of a transient southerly flow with reversals on 
a time scale as short as a day. 

Based on the present data and the previous works cited, the 
southward flow appeared to be of a transient nature and southerly 
origin. 

(5) The free-fall, STD method provided the S}moptic measurement 
of the velocity and density fields necessary for a study of the 
baroclinic structure of the Florida Current. Due to the tidal aliasing 
problem, time averaged mean values must be used. 

There are numerous other techniques that could have been used to 
process and analyze the free-fall, STD data. Some examples are: use 

of some other order of polynomial besides cubic; use of the spline 
technique for curve fitting rather than least squares technique; 



fitting a polynomial to the raw (i.e^ the actual values of transport 
measured) data to obtain a mean transport curve; and direct differen- 
tiation of the mean, or individual, transport curves rather than of 
the polynomial representative of the curves. Additionally, the same 
curve fitting technique used with the observed transport curve could 
be applied to the geos trophic transport curve. This would make the 
comparisons between the computed and observed quantities more uniform 
Because of the uncertainties in determining the **best” technique 
more emphasis in the future should be given to statistical and error 
analyses. Despite the lack of these analyses, the results of this 
thesis gave a realistic nine day mean description of the baroclinic 



structure of the Florida Current. 
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The average crossing time for a transect was about 8 hours. 

Because of the distorting influence of the periodic tidal forces, it 
was necessary to work with mean values averaged over the nine days 
rather than daily values. As shown in Table II, the station times 
were varied as practicable by starting a transect at Miami and Bimini 
on alternate days in an effort to reduce the tidal aliasing. 

Because of the short record length, a tidal analysis was not 
possible. The following method was used to determine the tidal aliasing 
effect . 

For the diurnal tide, consider 



where 

cv(3)is the normalized error in the mean, computed over N samples 
due to a diurnal (semidiurnal) constituent ; 




N 



and for the semidiurnal tide 




g _ 2 tt 

27 . 3 ( 24 ) 



lunar declination frequency (radians /hour ) ; 



6 - — 

~ 365 ( 24 ) 



solar declination frequency (radians /hour ) ; 



a = ^ 
25 



lunar tidal frequency (radians /hour ) ; 




solar tidal frequency (radians /hour ) ; 



N is the number of observations for a particular station; 



tj is the time of each observation for a particular station; 



and, 0 is the phase angle. 




i 



The phase angle 0 was determined from the difference between the value 
of equilibrium argument (Vq + u) of a constituent when t = 0 and the 
epoch of a constituent (k). The values of (Vq + u) were determined 
from Schureman (1958) and values of < from Smith, Zetler and Broida 
(1969). 



Using the four constituents 0 ^^ M 2 and 82 ? the velocity error, 

Av, of the mean observed velocity due to the diurnal (lur.ar, solar) 
and semidiurnal (lunar, solar) tides was computed: 

^^Diurnal(lunar ,solar) “ ^ ’ 

% 



and 



^^Semidiurna 1( lunar , solar ) ^ ’ 



"0 



where 

are the factors described previously, 

V-p is the current speed of a tidal constituent determined from 
Smith, Zetler and Broida (1969), 

and, Vq is the ’’average" current speed of the water column estimated 
by taking % of the mean observed surface velocity. 

The values of the parameters used in the calculations and the 
resultant values of % velocity error (AV) are shown in the tables below 
Tidal constituent 0]^ M 2 S 2 

Epoch of constituent (<) 24^ 12^ 115^ 309^ 

Velocity of constituent 
(Vp) cm/sec 



5.7 



5.6 



3.4 



2.5 



57 



Station Number 


1 


2 


3 


4 


5 


6 


7 


Velocity of water 
column (V-p) cm/sec 


33 


68 


92 


98 


96 


89 


80 


Normalized error (cv) , 
diurnal, lunar 


-.15 


-.27 


-.30 


-.27 


-.37 


-.41 


-.40 


Normalized error (cv), 
diurnal, solar 


.51 


.59 


.67 


.72 


.80 


.84 


00 

• 


Normalized error (P), 
semidiurnal, lunar 


.03 


.03 


.07 


.06 


.08 


.06 


.12 


Normalized error (3), 
semidiurnal, solar 


-.31 


-.17 


-.08 


- .06 


.17 


.32 


.40 


7o Velocity error (AV) 


4.3 


2.3 


2.4 


2.7 


3.3 


3.9 


4.9 


Station Number 


8 


9 


10 


11 


12 


13 




Velocity of water 
column (Vp) cm/sec 


63 


54 


51 


45 


42 


41 




Normalized error (Q^), 
diurnal, lunar 


-.36 


-.21 


-.22 


-.17 


-.17 


-.04 




Normalized error (cv) , 
diurnal, solar 


* 

00 


.70 


.49 


.35 


.33 


.26 




Normalized error (3), 
semidiurna 1 , lunar 


CM 

• 


,02 


.00 


.00 


-.01 


-.03 




Normalized error (3), 
semidiurna 1 , solar 


.24 


.21 


-.18 


-.27 


-.32 


-.23 




7o Velocity error (Av) 


5.2 


6.6 


2.1 


0.9 


0.1 


1,9 





Hence, the error in menn observed velocity values due to tidal 
biasing is about 37o. Also, as anticipated, errors due to diurnal tides 
predominate those due to semidiurnal tides and the error is greatest 



in mid-channel. 



m 
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